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We calculate the excitation spectra for the one-d 
Heisenberg-Ising antiferromagnets by expansions around the 
Ising limit. For 5* = 1/2, the calculated expansion coefficients 
for the spinon-spectra agree term by term with the solution of 
Johnson and McCoy. For S = 1, the solitons become gapless 
before the Heisenberg limit is reached, signalling a transition 
to the Haldane phase. By applying a staggered field we calcu- 
late the one-magnon spectra for the S = 1 Heisenberg chain. 
For 5 = 3/2 the quantum renormalization of the spin- wave 
spectra is calculated to be approximately 1.16. 



The study of spin-chains has produced many exciting 
ideas in many-body physics. The combination of exact 
solutions, field theories and numerical studies have led 
to a rather comprehensive understanding of the physical 
properties of these systems . One of the most remark- 
able physics to come out of one-dimensinal antiferromag- 
nets is the difference between the excitation spectra of in- 
teger and half-integer spin-chains. The S —1/2 chain is 
integrable and a lot is rigorously known about its excita- 
tion spectra, des Cloiseaux and Pearson |^ were the first 
to obtain expressions for the spectrum of the Heisenberg 
model. Later Johnson and McCoy obtained expres- 
sions for the low lying excitations of the Heisenberg-Ising 
model. 

Fadeev and Takhtajan Q| argued that the elementary 
excitations for the spin-half case are spinons and have the 
quantum numbers of a neutral fermion. Using semiclas- 
sical arguments, Haldane [|j showed that the topological 
excitations for integer spin-chains had integer spin and 
those for half-integer spin-chains had half-integer spin. 
He further argued that for integer spin these solitons 
become gapless before the Heisenberg limit is reached. 
Condensation of these solitons leads to a novel ground 
state, with a gap in the magnon excitation spectra. For 
half integer spins the spectra remains similar to the spin- 
half case, becoming gapless and linear at the Heisenberg 
point. By now there are numerous numerical |^ confir- 
mations of these remarkable facts. 

In this paper we show that the spectra for the topo- 
logical excitations such as spinons and solitons as well as 
the small spin-deviation excitations or magnons can be 
calculated systematically by a straightforward Raleigh- 
Schrodinger perturbation theory. This is simply a conse- 
quence of the fact that in one-dimension a domain-wall 
between two ground states is a localized object and has 



a well defined dispersion. While the perturbation expan- 
sion gives essentially the complete answer for systems 
with significant Ising anisotropy, it is also very accurate 
for Heisenberg systems. For S = 1/2, our calculated co- 
efficients for the spinon spectra agree term by term with 
a direct expansion of the closed form results of John- 
son and McCoy We also calculate the soliton ex- 
citation spectra for the S = \ chain and find that it 
becomes gapless before the Heisenberg model is reached. 
As shown by Haldane these solitons play the same 
role as the domain-walls in the 1 + 1-dimensional classi- 
cal Ising model, where vanishing free-energy cost for the 
domain-walls leads to a second order phase transitions. 

In addition to the topological excitations we also calcu- 
late the spectrum for magnons which are the quanta for 
small spin-deviations around a given state. For S = 1/2 
chain, a magnon is not a sharply defined state even in the 
Ising limit as its energy is degenerate with that of two 
spinons. Thus it can break apart into two spinons. For 
S > 1, magnons are well defined elementary excitations 
near the ising limit. For = 1, we add a staggered field 
term to the Hamiltonian, with a coefficient which goes 
to zero at the Heisenberg point. The presence of this 
term avoids the condensation of solitons with non-zero 
Ising anisotropy, and allows us to calculate the magnon 
excitation spectra in the Heisenberg limit. Even simple 
addition of available terms in the series leads to a spec- 
trum which agrees with previous numerical estimates to 
within a few percent. 

For = 3/2 we are not aware of any previous numer- 
ical calculation of the spin-wave excitation spectra. Our 
calculations provide extremely accurate estimates of the 
spectra for systems with significant Ising anisotropy. We 
find that the series also show good convergence upto the 
Heisenberg limit especially away from the gapless point 
near q = Q, where they are consistent with a linear spec- 
tra but converge more slowly. We estimate the quan- 
tum renormalization of the spin- wave spectrum from the 
maximum spin- wave energy as was done in recent exper- 
iments on CsVCls . We obtain this renormalization to 
be approximately 1.16. This is larger than the order 1/ S 
spin- wave estimate of 1.12 ||] but still smaller than the 
central estimates of the experiments which give 1.26 0. 
However, the experimental uncertainties are ±0.2 due to 
large uncertainties in the value of J. 

In the Ising limit all these excitations are localized, and 
those localized at different points in space are degener- 
ate. The magnons consist of a single spin-deviation in a 
Neel state, whereas a soliton or a spinon is a domain wall 
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between two degenerate ground states. A simple way 
to think of the topological excitations is in terms of the 
ground state of a large odd-site ring. In the Ising limit, 
an Ai"-site ring has at least N degenerate ground states, 
corresponding to the different locations where the neigh- 
boring spins are not antiparallel. This broken bond can 
be thought of as the domain-wall. Since, these ground- 
state configurations differ locally from the ground state 
of an even site ring only at the wall, ie the energy-density 
with respect to the ground-state is non-zero only on one 
bond, we can regard the excitation as localized there. As 
the X — y coupling is turned on these degenerate ground 
states evolve into a band of states with well-defined en- 
ergy momentum relation. Thinking in terms of an odd- 
site ring also explains why these domain-wall excitations 
have half-integer spin for spin-half chains and integer spin 
for integral spin-chains. In the Ising limit a spin-chain 
with S > 1/2 has other ground-state degeneracies, for 
example the spin right next to the broken bond can be 
in any one of the states. Most of the degeneracies are 
lifted with non-zero x — y coupling and only those related 
to symmetry survive. 

The expansion coefficients for the excitation spectra 
are obtained by a degenerate perturbation theory which 
has recently been formulated in terms of a cluster expan- 
sion by Gelfand |^ . The magnon calculation is similar to 
that in higher dimension jlOj , but the solitons are special 
to one-d as discussed in the previous paragraph. For the 
purpose of our calculation we only need to focus on the 
vicinity of the excitation. The boundaries can be con- 
sidered to be infinitely far away and do not matter. The 
degenerate perturbation theory, allows us to construct an 
effective hamiltonian in the space of these localized Ising 
states, which can then be diagonalized straightforwardly 
by Fourier transformation. The recursive procedure for 
obtaining the effective Hamiltonian has been discussed 
by Gelfand §. 

We consider the Heisenberg-Ising Hamiltonian: 

n = jY^ sfs!^, + \{sfsf^, + 5f (1) 

i 

For the 5 = 1/2 case it is most convenient to present the 
square of the spinon dispersion, which we calculate via 
the cluster expansion to order A^^ 

lQe{qf =4 + 20A2-6A4 + |A6 + |A«-y|gAi2 

-l-cos2g(16A + 4A3 - \\'^ - ^^A^ + ^A") 

+ 0(Ai3) (2) 

A closed form expression for the spinon dispersion can 
be read off from the work of Johnson and McCoy |^ : 

e[q) = (1 - A2)1/2-(1 - fc2 sin2 qf/^ (3) 

TT 

Here K ( and similarly if' ) is a complete elliptic integral 
of the first kind with modulus k. The latter is determined 
as a function of A from the relation 



K'jk) ^ cosh-^(l/A) 

K{k) TT ^ ' 

We have verified that to the order calculated our expan- 
sion coefficients agree term by term with a direct expan- 
sion of the above expressions In general, this expan- 
sion shows excellent convergence. By simply adding the 
terms in Eq. 2, one obtains the dispersion for the Hciscn- 
berg model as shown in Figure 1 and compared with the 
exact Heisenberg spectrum. Except very near q = 7r/2, 
where the spectrum becomes gapless at the Heisenberg 
point, the convergence is excellent. It is also evident that 
the maximum of the spectra at g = is much less sensi- 
tive to small anisotropy than the spectrum near q = 7r/2. 

For 5 = 1, the soliton states can have 5^ = 0, ±1. 
While all three are degenerate in the Ising limit, the 5^ = 
soliton becomes the lowest energy one for A 7^ 0. We 
have calculated the full spectrum for the solitons. The 
minimum energy soliton corresponds to g = tt, unlike the 
spin-half case where the minimum energy is at g = 7r/2. 
The mass gap is given by: 

m = 2 - 2A - f A2 + |a3 - 1.207407A* + 1.971667A5 
- 2.723225A6 + 3.294886A^ - 3.400478A8 

-f 2.240766A9 + 1.571335A10 + ... (5) 

Since we expect a 2D Ising critical point, with v = 1, 
where the soliton becomes massless, we use simple Fade 
approximants to obtain the location of this critical point. 
We estimate Ac = 0.842±0.002, which compares well with 
previous studies of this phase transition using ground 
state properties p^ . 

In order to calculate the magnon spectra in the Heisen- 
berg limit we consider the Hamiltonian: 

+ (l-A)E.(-iy^f- (6) 

The addition of the last term ensures that for A ^ 1 , only 
one of the two Neel states is the ground state of the sys- 
tem and the other Neel state ( and the solitons) have infi- 
nite excitation energies. Exactly at the Heisenberg point 
( A = 1 ) the staggered field term goes to zero and one 
recovers the Heisenberg hamiltonian. Because we have 
introduced the staggered field and doubled the unit cell, 
the magnon spectra is also doubled. Since there is no long 
range order in the Heisenberg limit, only the calculated 
spectra between wavevectors (7r/2, tt) will survive and the 
spectral weights for those between (0, 7r/2) will vanish as 
A 1. The effective Hamiltonian, whose fourier trans- 
form gives the excitation spectra is given complete to 
order A^" in table 1. The convergence of the series for 
A = 1 is excellent. In Figure 2 we show simply the sum 
of terms from 8, 9 and 10 term series. Fade approximants 
can be used to improve convergence slightly. They lead 
to estimates for the energy gap of 0.42±0.01J. The entire 
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spectrum is within a few percent of those calculated ear- 
lier by exact numerical diagonalization of finite systems 

For S — 3/2 we have also calculated the effective 
hamiltonian for the magnon excitation spectra. It is pre- 
sented in table 2. In Figure 3 we show the estimated 
dispersion for A = 1 obtained by simply adding terms 
upto 4th, 6th and 8th order. Except very near q — 0, 
where the dispersion should become linear and gapless 
at the Heisenberg point, we find that this simple addi- 
tion converges rapidly. By assuming a linear dispersion 
near q — 0, and using series extrapolation methods we 
can directly estimate the spin-wave velocity. However, 
if we assume that the spin-wave spectrum is uniformly 
renormalized with respect to the classical spectrum, the 
quantum renormalization of the spin- wave spectrum can 
be estimated from the value of the excitation energy at 
q — 7r/2. The latter measure, involving the maximum 
spin-wave energy, is insensitive to anisotropy, whereas 
the linear spectrum at small q is clearly very sensitive 
to anisotropy. Using the latter measure, we estimate the 
quantum renormalization of the spectrum to be approxi- 
mately 1.16. This is somewhat larger than the order l/S 
spin- wave estimate of 1.12 Recently, this quantum 
renormalization of the spin-wave spectra has also been 
determined experimentally in the material CsVCls from 
the energy of the zone-boundary magnons Q . They ob- 
tain a spin-wave renormalization of 1.26 ± 0.2, the large 
uncertainties resulting from the lack of accurate determi- 
nation of the exchange constant J. 

In conclusion, we have shown that topological excita- 
tions such as spinous and solitons as well as magnons in 
one-dimensional Heisenberg-Ising antiferromagnets can 
be studied by a straightforward application of Raleigh- 
Schrodinger perturbation theory. This expansion gives 
nearly the complete answer for systems with signifi- 
cant Ising anisotropy, but is also highly accurate for the 
Heisenberg model. We have used this method here to 
verify the spinon spectrum for S = 1/2, to calculate the 
soliton spectrum for 5 = 1 and show that it becomes 
gapless at the Ising critical point, as well as to calculate 
the magnon spectra for S — 1 and S = 3/2 chains. The 
application of a staggered field, whose coefficient goes 
to zero in the Heisenberg limit allows one to study the 
magnon spectra for the Haldane-gap Heisenberg systems 
as well. 
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FIG. 1. The spinon excitation spectra for the spin-half 
Heisenberg model obtained by adding terms upto order 10 
(crosses), 11 (circles) and 12 (squares) in Eq. 2. The solid 
line is the exact result 

FIG. 2. The magnon spectra for the spin-1 Heisenberg 
model obtained by adding terms upto order 8 (squares), 9 
(circles) and 10 (crosses) in table 1. 

FIG. 3. The spin-wave spectra for the spin-3/2 Heisenberg 
model obtained by adding terms upto order 4 (crosses), 6 
(circles) and 8 (squares) in table 2. 

TABLE I. Effective Hamiltonian for the 5 = 1 Heisen- 
berg-Ising model (Eq. 6). To get the magnon dispersion at 
wavevector q one needs to sum the series over r with a factor 
cos(gr). 

r Series 

3 - A - 0.2A^ - 0.18A'' - 0.0530357A'' + 0.0088342A'' 

-l-0.0109158A^ + 0.0015458A^ + 0.0021357A* 

-f 0.0080557A^ + 0.0105825A1'' 
2 -iA^ - iA^ - 0.1094048A'* - 0.0644654A^ - 0.0528568A'^ 

-0.0381666A^ - 0.0176628A* - 0.0020513A^ 

-F0.0025792A" 
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4 -O-OigSTSA"* - 0.0245625A^ - 0.0212035A'^ - 0.0177130A'^ 
-0.0153450A^ - 0.0132463A'^ - 0.0114761A^° 

6 -0.0016698A'' - 0.0035734A'^ - 0.0046995A® 
-0.0052158A® - 0.0054470A^° 

8 -0.0001682A* - 0.0004937A^ - 0.0008462A" 

10 -0.0000195A^° 



TABLE II. Effective Hamiltonian for 5 = 3/2 Heisenberg 
Ising model. To get the magnon dispersion at wavovcctor q 
one needs to sum the series over r with a factor cos{qr). 

r Scries 

3 - 0.825A^ - 0.0529051A* + 0.0195540A'' + 0.0186877A'* 

2 -1.125A2 - 0.2552009A® - 0.0900226A® - 0.0174543A® 

4 -0.0980859A'' - 0.0691161A^ - 0.0454369A* 

6 -0.0190201A'' - 0.0229853A* 

8 -0.0043107A® 
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